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ABSTRACT 

A family of cusped potential-density pairs is presented for modelling galaxies and dark haloes. 
The density profile is cusped like p ~ r''^- at small radii. The distribution function is simple 
and takes the form / oc j^p-^E^^P+^y^ (where E is the binding energy and L is the angular 
momentum). The models all possess the remarkable property that the virial theorem holds lo- 
cally, from which they earn their name as the hypervirial family. Famously, this property was 
first discovered by Eddington to hold for the Plummer model in 1916. In fact, the seductive 
properties of the Plummer model extend to the whole hypervirial family, including the mem- 



bers possessing the cosmologically important cusps of p 



r ' orp 



r or p 



r-^l\ The 



intrinsic and projected properties of the family of models are discussed in some detail. 

Key words: stellar dynamics - celestial mechanics - galaxies: elliptical and lenticular 
galaxies: kinematics and dynamics - galaxies: structure. 



1 INTRODUCTION 

The allure of spherical models is their simplicity. Certainly, there 
are few exactly spherical galaxies or bulges. However, spheri- 
cal models are still useful as representations of galactic nuclei, 
where the kinematics are usually dominated by a central cusp (see 
e.g., Dehnen 1993; Tremaine et al. 1994). Almost nothing is firmly 
known concerning the shapes of the dark matter haloes around disc 
galaxies. So, it is very reasonable to use spherical models to inves- 
tigate, for example, the kine matics of distant halo stars and satellite 
galaxies of the Milky Way ('Whitell 19851: iKulessa & Lvnden-Belll 
[l992; Wilkinson & Evans 1999). Globular clusters with long relax- 
ation times, such as (jj Centauri, roundish dwarf spheroidal galaxies 
such as Draco ("e.e.. iWilkinson et al.,.200Z) . and clusters of galaxies 
are also possib l e appl ications for spherical models. 

IPlummeJ il91ll) set the standard when he introduced the 
model with potential i/r and mass density p as a fit to the outer parts 
of globular clusters. 
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The real beauty of t his model lies in two facts discovered - not 
by Plummer - but bv lEddinetoJ Jl9ld) . First, the isotropic Plum- 
mer model has a very simple distribution function. In fact, Edding- 
ton showed that the isotropic distribution functio n of any sp herical 
model is a vailable immediately as a quadrature teddingtonll 1 9 1 d 
iBinnev (^Tremaine. 1987. p. 237), 



f(E) = 
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Here, / is the distribution function, which depends on the specific 
binding energy E of the stars alone. Substituting the particular case 
of the Plummer model, Eddington found the very simple answer 
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Second, lEddingtonlll91^ also pointed out that the isotropic Plum- 
mer model has the remarkable property that it obeys the virial 
theorem locally (of course, any self-gravitating system must obey 
it globally). At every spot, the kinetic energy in each element 
(r = p(v^)/2) is exactly one-half of the magnitude of the local 
potential energy (W = -pi/'/2). Therefore, the virial theorem holds 
(2T -I- W = 0) at each and every spot! Such models we call hyper- 
virial. 

In this paper, we introduce a new family of cusped galaxy 
models. Each member of the family has a distribution function that 
rivals equation in its simplicity. Even more remarkably, each 
member of the family satisfies the virial theorem locally. The Plum- 
mer model is simply the limiting case when there is no cusp. 



2 INTRINSIC PROPERTIES 

The potential tf/ and density p of the family of models are 
GM (p+l)M a'' 



(a'' + rPyp' 
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where p is a parameter that is a positive real number. The model 
is the limiting c ase of the ' generaliz ed isochronous' models first 
written down bv I VeltmannI il979ald) . who also worked out some 
of their integrated properties. He, however, presented the potential- 
density pair as a useful fitting formula and did not study its 
dynamical basis. 
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The case p = 2 is recognized as the Plumm er model, while th e 
case p = I corresponds to a model introduced bv lHemauistl<199d) . 
Generally, the density distribution behaves like r"<^"''> at small radii 
and falls off like r"*''^^' at large radii. Note that, if p > 2, there is a 
hole at the centre |p(0) = 0] and the density increases outwards near 
the centre. The model with p = cx) is a shell of mass M and radius 
a. Models with p > 2 are therefore not astrophysically realistic. 
Unless otherwise noted, we only consider models for which the 
parameter p is restricted to lie in the range < p < 2. Henceforth, 
we use units in which G = M = a = 1. 

The cumulative mass within the sphere of radius of r is 



1 



(5) 



V = 47t f p/Ar = -r^^ = 

Jo dr (1 + ri'Y+^'P (1 + r-py+^'P 

and so the half-mass radius ri/2 is 

= [2"«"^" - 1]-"". (6) 



The circular velocity curve Vdrc is 
2 dif/ M, rP 



dr 
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(7) 



As r — » 0, the circular velocity tends to zero. As r oo, the cir- 
cular velocity becomes asymptotically Keplerian. The velocity dis- 
persion is determined by solving the Jeans equation 



-:r[p(v,)) + 2l3 = — 

p dr ^ ' r or 



(8) 



where f} is the Binney's anisotropy parameter I - p = (v^>/(2(v^», 
and (v^> and (v^) are the squares of the radial and tangential ve- 
locity dispersions. In general, for a given potential-density pair, if 
the behaviour of /3 is assumed, then the Jeans equation Jsj can be 
solved using an integrating factor. In particular, if /3 is a constant, 
then the radial velocity dispersion becomes 
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where 2Fi(a,b;c; x) is the Gaussian hypergeometric function. In 
general, (v^ > is finite everywhere if 2yS < (2-p), while it diverges at 
the centre if 2yS > (2 - p). In fact, we find that the models with 2fi > 
(2-p) are unphysical as the distribution function is not everywhere 
non-negative ( Appendix IXI see also .An & Evans..2005.') . 



2.1 Hypervirial Models 

It is of course awkward to work with hypergeometric functions, 
and so it is natural to look for simplifications. Fortunately, a won- 
derful simplification exists. If 2/3 = (2 - p) in equation J9j, then 
the hypergeometric function becomes the constant unity! So, the 
square of the velocity dispersion becomes linearly proportional to 
the potential everywhere, i.e.. 



2p + 2(l + rPy/P 2(p+iy 

pif/ 
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Furthermore, we deduce the unusual result that the total velocity 
dispersion, (vj) + (v|) = tf//2, is independent of p. This leads us 



to suspect that the distribution functions of these models must also 
be very simple, though they must be dependent on the angular mo- 
mentum L and the binding energy E since the velocity dispersion 
tensor is in general anisotropic. It is straightforward to show that if 
the angular momentum dependence of the distribution function is in 
the form of a power law, i.e., f(E, L) = L^^^fE(E), then the veloc- 
ity dispersion anisotropy is everywhere const ant and t he Binney's 
param eter becomes the constant j8 (see e.g., iBinne v & Tremain^ 
Il987ft . Motivated by this, we expect the distribution functions of 
our models will also have the form of /(£, L) = LP^^ fE(E), where 
fE(E) is a function of energy that remains to be found. 

A little more work shows that fE(E) is itself a power law. This 
follows because the density can be decomposed in terms of a prod- 
uct of powers of the potential and radius 
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and consequently simple distri bution functions exist of the form 
(see e.g., eq. B7 of lEvanJl994l) 



f{E,L) = CLP-^E'^p^^^''-. 



(12) 



Note that f{E, L = 0) diverges for p < 2 which is in accor- 
dance with the presence of the density cusp at the centre. Here, 
the constant C can be determined by the normalization condition 
p = //d'v and is 

1 r(2p + 3) 



C 



(13) 



2;'/2+i(27T)5/2 r(p/2)r(3p/2 + 3/2) 

where r(.i:) is the Gamma function. The velocity dispersion ten- 
sor is found from the second velocity moments of the distribution 
function. 
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(15) 



When p = 2 (the Plummer model), the velocity dispersion is 
isotropic and the distribution function no longer depends on the 
angular momentum. As p — > 0, the density of the models become 
increasingly cusped and the velocity distribution becomes increas- 
ingly dominated by radial orbits. 

Note that the kinetic energy T in any element at each point is 
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where W = -(ptf/)/2 is the local contribution to the potential en- 
ergy. In other words, we have established that there exists a virial 
relation 2T + W = that holds at every spot. All these models are 
therefore h ypervi rial. This generalizes the remarkable result that 
lEddingtoi ? ( 1916) originally found for the Plummer model. 

In fact, Baes & Deionghe ( 2002) have already noted that the 
Hernquist model with /3 = 1/2 can be be constructed from the dis- 
tribution function 



f(E,L): 



3 E- 



(17) 
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This is a particular case of equation <12t when p - 1. In 
fact, all the models described by equation <12t c onstitute spe- 
cial cases of the 'generalized polytropes' (see e.g.. iHenorll 197.31 : 
iBarnes. Goodman & Hut|[l986l) with their potential expressible as 
elementary functions of the radial distance. As we show in the next 
section, they are also the only generalized polytropes with a finite 
mass and an infinite extent. 



2.2 Power-law distribution functions 

Now, let us ask the question: are there any more spherical models 
with simple distribution functions that are just a power of energy 
multiplied by a power of angular momentum? 

Suppose that the distribution function of a spherically sym- 
metric system is indeed given by the ansatz 

f(E,L) = CL-^'^E"-^'^ (18) 
where yS < 1 and n > 1/2. Then, the corresponding density becomes 



r(n-/3+ 1) r^P ■ 
The velocity dispersions may also be found as 

r(l -/3)r(;2- 1/2) 
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However, the total kinetic energy T^tot 

r„ = i/p«v?) + <v4))dV = i/dVp^^ 

3-2/3 



(20) 
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where Wtot is the total potential energy. In other words, the sum 

W,„, + 2T,„, is 



1 - 



2(3 - 2/3) 
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(22) 



suggesting that the global virial theorem for a steady-state system 
is satisfied only if ;i + 3j8 = 5 [or, in other words, if p = 1 - {p/2), 
then n = 2 + (3/2)p and n - 3/2 = (3p+ l)/2]. In this case, the total 
potential energy can be explicitly evaluated as 



-2r,o, = -27t I pi//r dr 
Jo 
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p = l/2 

Hernquist (p — 1) 
p=3/2 

Plummer (p = 2} 




Figure 1. The surface density profile of the hypervirial models with p 
1/2, 1 (Hernquist), 3/2 and 2 (Plummer). 



where p = 2(1-/3) = 2(n- 2)/3. So, our hypervirial models are the 
only spherically symmetric steady-state systems with finite mass 
which have distribution functions as simple as equation <18> . In 
particular, the density profile of the Plummer model is the only 
spherically symmetric self-gravitating polytrope of finite mass with 
an infinite extent that can be thermally (pressure) supported. 

There are other known spherical models which do have distri- 
bution functi ons as simple as equation <18> - namely, the power- 
law spheres jEvanslll994l) . However, these models have infinite 
mass and so do not satisfy the global virial theorem (unless bound- 
ary terms are added). Similarly, stellar dynamical polytropes also 
exist, but have di stribution function s / oc (£ - Eq)'^, where Eq and 
A' are constants jEddingtorjfigT^ . Only for the Plummer model 
does £0 vanish and so the model is of infinite extent. 



3 PROJECTED QUANTITIES 

The beauty of the hypervirial models lies in their simple distribu- 
tions of velocities. In contrast, the projected quantities are generally 
more awkward, typically reducing to elementary functions only in 
the cases of Plummer (p = 2) and Hernquist (p = 1).' 
The surface density is 



J'""' prdr p -t- 1 „ , r"'' CO 
« V72^ 27t Jo (R" + 



;P+'0de 



cos''6')2+i/''' 



(24) 



If the mass-to-light ratio is constant, then this is proportional to 
the surface brightness. Fig. Q shows the surface brightness profile 
for a number of the hypervirial models. For the t wo special cases , 
expressions usin g elementary functions are known <Plummeill9Tll 
lHernauisJl990l) 



1 



7t(l +ij2)2 



Plummer (p = 2), 



' Although we do not derive the explicit forms, the quantities for the case 
of p = 1/2 in general reduces to the expression involving elliptic integrals. 
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p = i/a 

Hernquist (p = l) 
p=3/2 

Plummer (p-2) 




Figure 2. The surface brightness of the hypervirial models with p = 1/2,1 
(Hernquist), 3/2 and 2 (Plummer) plotted against (R/Re)'''* where R/ife is 
the projecte d dist ance in units of the effective radius of each model. In this 
diagram, the lde V aucouleurs ( 1948) profile would be a straight fine with a 
slope of -3.331. The inset shows the residuals with respect to the de Vau- 
couleurs profile. 



^ [(2 + R-) X(R) - 3] Hernquist (p = 1), 



2n (l-R^) 

where, the Hernquist function X(R) is defined as 

2 / 131-i?\ arcsech arcsec 



X(R)- 



1+R 



2 2 l+Rj 



Despite the appearance of the formal singularity al R = I for the 
p = I case, Eh is in fact regular everywhere for R > O.ln particular, 
it is continuous [Xh(1) = (27t)~'(4/15)] and differentiable a.lR= 1. 
The central surface density is finite if p > 1, 



E(0) = 2 r pdr: 
Jo 



r(i-i/p)r(i + 2/p) 
27tr(i + i/p) 



(25) 



but it is divergent otherwise. In particular, if < p < 1, the central 
surface density behaves like 



1 (p + i)r(i/2-w2) 



oc R-'^^-"^ (R 0), 



(26) 



Rl~P 47tl/2r(l -p/2) 

and it diverges logarithmically if p = 1 faemauisll990h . At large 
radii, the asymptotic behaviour of the surface density is given by 

1 r(i + p/2) 



27ti/2r(l/2 + p/2) 



oc R- 



(R 00). 



(27) 



If the mass-to-light ratio is constant, then the cumulative 
brightness is proportional to the mass within the cylinder of radius 
of R, namely 



Jo Jo (R" + cosi'ey+'i" 



(28) 



For the Plummer and Hernquist models, the result reduces to 
R'- 



WlP 



Plummer (p = 2), 




100 



R 



Figure 3. The line-of-sight velocity dispersion profiles for the hypervirial 
and circular orbit distribution functions corresponding to the models with 
p = 1/2, 1 (Hernquist), 3/2 and 2 (Plummer). 



l-R^ 



[XiR)-l] Hernquist (p = 1). 



Note that Ml(l) = 1/2 and M^(l) = 1/3. The efi'ective radius (or 
the half-light radius) R^ can, in general, be found by numerically 
solving 

Mr(R,) = 

The particular solutions are ~ 1 (exact), 1.18084, 1.81527 and 
11.0151 for p = 2, 3/2, 1 and 1/2, respectively. Fig. |5| shows the 
surface brightness normalized by E(i?e) as a function of (R/R^y^*. 
We find that Hernquist profile - and by extension the hypervirial 
models with p a; 1 - are all reasonably good approximations to the 
empirical de Vaucouleurs law ( 1948) of the surface brightness 
profiles between 0.25 < (R/Rc) ^ 10. The total mass is of course 



Mr(oo) : 



itm r ■ 

R-»Jo (1+^ 



cosede 



-p cospey+^'p 



C 
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cosSde = 1, 



as it should be! 

For any spherical model with arusotro py parameter ff, the 
line-of-sight velocity dispersion Vo^ is (e.g.. lBinnev & Tremaind 
1987) 



hypi 

^Rp- r 

47t Jo 



For the hypervirial models,/? = 1 - (p/2), and so 
cos"+20 -I- (p/2 - 1) cosP+^e 



(RP + COS''e)2+2/P 



A0. 



(29) 



(30) 
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As in the case of the surface density, this integral reduces to ele- 
mentary functions for the Plummer model 

3 1 , 37t 1 



^^'^^ 64(1 +i?2)5/2' 
and for the Hernquist model 



1 



64 (1 



- [(24/?* - 68/?-* + SIR- - 28) 



An { 12(1 -i?2)3 
+3(8/?* - 28i?* + - UR- + 4) X{R)^ . 



(31) 



This is continuous al R = \ with cr^(l) = [47tEH(l)]"' [tt - 
(326/10 5)]. Note that, the isotropic model was the main object of 
study of'Hernauist' ('1990*) and its line-of-sight velocity dispersion 
(eq. 41 of Hernquist 1990) dilfers from equation <31> . of course. 
If p > 1, then the central value of the velocity dispersion is 



O-(0) : 



)dr : 



r(i + 3/p)r(i + i/p) 



4r(2/p)r(2 + 2//7) 

which varies from 37t/64 for p = 2Xo 1/4 as p — > 1^ 
hand, if < p < 1 , the central velocity dispersion is 

2 1 (p+l)r(l/2-p/2) E 

~ 4 



Ri-P 167t'/2r(l - p/2) 
For the case p = I, using X(R) ~ ln{2/R) as R 



(R 0). 



(32) 

. On the other 
(33) 



' 0, we find that 



n\ R 2 



limcrr, = lim 



In other word, the line-of-sight velocity dispersion at the centre is 
always finite for < p < 2. In particular, cr^(O) = l/4if0<p<l 
while it can be determined from equation <32> for 1 ^ p < 2. At 
large radii, the fall-off becomes Keplerian 



2 1 (p + 2)r(3/2 + p/2)^ , 
cr ~ — ^ oc R 



{R oo). 



(34) 



Rsr{3 + p/2)r(i + p/2) 

It is instructive to compare this with the purely circular orbit 
model, for which the line-of-sight velocity dispersion is given by 



Jr 

p+1 



2 r2 - R2 



Rip- 



47t 

which becomes 

EpCTp 



Jrrc/l 
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157t 
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128 il+R^f'^-' 



128 (1 +i?2)3/2 



1 



nR- 



- [(24R' 



(35) 



(p = 2), 



154) 



47t [ 12(1 -i?2)4 

+3(8/?* - 36R^ + 63R* - 40R^ + 40) X(R)\] (p = 1). 

Since we only observe radially directed motion when looking at the 
very centre, the line-of-sight velocity dispersion must vanish there 
for the extreme tangentially anisotropic model. On the other hand, 
at large radii, we have 



1 r(5/2 + p/2)r(3/2 + p/2) 



oc R- 



(R oo). 



(36) 



R 2r(3 + p/2)r(l + p/2) 

By comparing this to equation <34t . we find that the line-of-sight 
velocity dispersion for the purely circular orbit model is larger than 
that for the hypervirial model at large radii. This is an expected 




hyperviria 



Figure 4. The line-of-sight velocity dispersion profiles for the isotropic, 
hypervirial and circular orbit distribution functions corresponding to the 
Hernquist model (c.f., fig. 1 of Hernquist 1990). The inset is the same but 
as a function of log R, which shows the behaviour near the centre better 



result, since we observe tangential motion preferentially as the lines 
of sight moves to the outskirts of the system. In Fig. we plot 
the line-of-sight velocity dispersion for the hypervirial and circular 
orbit model for a few representative values of p. In Fig.|4| the line- 
of-sight velocity dispersions for the Hernquist model are shown. 
The expression for the velocity dispersi on of t he isotropic model is 
provided by equation (41) of Hern guisl il990l) . 



4 HYPERVIRIALITY 

Next, let us ask the question: are there any more spherical models 
with the property of hyperviriality? Mathematically, such a model 
has to satisfy 



<v?> + <v4) 



(37) 



Since 1 - /? = (v^>/(2(vj», if f} is constant, equation <37t implies 
that 



(38) 



2(3 - 2/3) 

On the other hand, the Jeans equation Jsj can be written 

'-^A-T^-V^'^^-T- (39) 

ar \p or r j ar 
From equations <38> and <39> . by eliminating (vj), we obtain 



dp 
dr 



i^5-Ap)'^-2p 
dlnr 



(40) 



Now, Poisson's equation for a spherically symmetric system reads 



1 d 



,d^ 



(41) 



where the negative sign is due to our choice of the sign for the po- 
tential. By eliminating p from equations <40t and <4H . we can get 
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a third-order non-linear differential equation for ip. However, using 
the relation di/r/dr = if/(dlm///dr), it is possible to reduce the equa- 
tion into a second-order differential equation for;^^ = (dlni/'/dr). 



— + 2^ 

dr^ dr 



I +13 



■ (1 - 2l3)x 



-2 



(1 -2y3)4 + (4 - 5yS)— + 2(1 



(42) 



or additionally changing the independent variable into f = In r, 

^1 - (1 - ipm + 1) J - 2(1 -pm + m + = o (43) 

dr dt 

where ^ = (dlni/f/dlnr) = (dlni/r/dO. 

Since equation 1431 does not involve the ind ependent va riable, 
its order may be reduced by the substitution (e.g. llncell944h 



d^ d-if_d^_d^d^ 
dF' df2" ~ dF ~ d^ df 

so that one finally arrives at 

2pm + i)f - 2(1 -pm + m + = o. (44) 

Since this is a first-order equation, it is always possible (at least 
formally) to find an integrating factor, which in this case is 

1 



Then, the equation can be written in the exact form 
^ In |[f - 2(1 -pm + l)]"'""' + «^ + 1)]| = 0. 
If yS < 1 , we obtain the solution 

u - 2(1 -pm + 1)]'""'" u + «f + D] = c 



(45) 



where C is constant. If we restrict attention to systems of finite 
mass, then the potential is asymptotically ~ r"' and consequently 
we can set C = from the boundary condition at r = oo (^ = - 1 
and ^ = 0). Then, we find two possible solutions, namely, 

^ = 2(l-y3)^(^-H), ^ = -^(^+1). 

Note that the second solution is independent of p. In fact, we find 
that it leads to vanishing density everywhere so that the second so- 
lution is unphysical. For the first solution, we obtain 



dt 



pdt 



where p = 2(1 - /?). Then, integrating gives the result that 



In 



f+1 



pit - to), e, ■■ 



1 



1 : 



(46) 



where the double sign in front of the exponential appears when the 
absolute value is removed. However, since f = (dlni/^/df), one can 
integrate equation H6i further to find i//, that is. 



In 



i/'o 
so that 



-It 



dt 



1 



-;'('-'o) 



■In 1 : 



py'('-'o)| 



(47) 



i-l/P 



(48) 



^ ^ L J [1 ± (r/ro)"]i/f 

Note that the absolute value in equation 147 > can be ignored since 
1^0 can be either negative or positive. The solution corresponding 
to the negative sign in equation <48> is unphysical, while the other 



solution recovers our hypervirial models. One integration constant 
(t(i = In ro) corresponds to the scalelength and the other integration 
constant (if/o) becomes the overall scaling factor. 

Our original differential equation before reduction was in third 
order, so we expect the general solution to contain three constants 
of integration. The hypervirial family is a two-parameter family, 
and so there exists an additional family of solutions, corresponding 
to the choice C in equation <45> and consequently having either 
infinite mass or finite extent. Unfortunately, the potential-density 
pair now seems not to be expressible in terms of elementary func- 
tions. 

Physically speaking, hyperviriality is related to stability of the 
model against evaporation. The escape speed is ^f2i//. In a hyper- 
virial model, the root mean square speed is ^fpjl. For any other 
kind of model, the root mean square speed may lie below the hy- 
pervirial value at some places, provided that there is suitable com- 
pensation at other places so that the global virial theorem is obeyed. 
In other words, at some spots, there will be more stars close to the 
escape speed and hence on the verge of escaping. If a few stars 
escape by accident or tidal perturbations, then the potential is low- 
ered and stars originally safe would be left with speeds above the 
escape speed. Hyperviriality therefore aids stability by minimizing 
the number of high-velocity stars. 

However, evaporation is not the only cause of instabil- 
ity. Dynamical effects, such as the radial orbit instability, may 
drive spherical stellar sys tems into triaxiality. In particular, 
IPalmer & Papaloizoii il987h argue that any distribution functions 
unbounded at zero angular momentum, as are all those given by 
equation <12t except for p = 2, are formally unstable, though the 
growth rate of any instability may be slow. Simulations of general- 
ized polytropes hav e already been carried out bv .Henon. (.1973^ and 
iBarnes et alJ il98d) . Such numerical work has tried to identify a 
critical ratio of kinetic energy in radial to tangential motion IT^/Tj 
above which the radial orbit instability occurs. This quantity is 2/p 
for the hypervirial models and so is 2 for the Hernquist model. 
The stability criterion 2T,./Tj < 2 has been suggested as a crude 
summ ary of a wide range of numerical experiments Bar nes et alJ 
Probably, only full-scale numerical simulations of the hy- 
pervirial models can find the exact point at which the radial orbit 
instability sets in, but it seems reasonable to expect the most highly 
cusped members (which are the most radially anisotropic) to be 
susceptible to the radial orbit instability. 



5 CONCLUSIONS 

This paper has provided a set of very simple distribution functions 
for a family of cusped spherical galaxy models. If the density is 
cusped like p ~ r''^^ at small radii, then there is a simple anisotropic 
distribution function, which behaves like / oc /^p-2£-(3p+i)/2 (^vhere 
E is the binding energy and L is the angular momentum). We call 
these models the hypervirial family. This is because every model 
obeys the virial theorem at each and every spot, and, of course glob- 
ally. 

The family includes the Hernquist model which possesses the 
cosmologically important r"' cusp at small radii and has a sim- 
ple distribution function / oc L"'£-. As its sole isotropic and un- 
cusped representative, the family includes the Plummer model with 
its familiar distribution function / oc There are also members 
which possess other density cusps such as p ~ r^"^!^ and p ~ r^^^~, 
which have been suggested as important on cosmogonic grounds 
cEvans & ColletL 1997; .Moore et al. 1998 ). 
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In a sense, the models in this paper really are the last word 
in simplicity - for we have proved that they are the only self- 
gravitating, spherically symmetric finite mass models with such 
simple distribution functions. 
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APPENDIX A: CONSTANT-/3 DISTRIBUTION FUNCTION 

Systems with a constant velocity anisotropy parameter always can 
be built from the distribution functions of the form 

f(E,L) = L-'-l'fdE) 

where fcCiS) is a functi on of the binding energy alone (e.g., 
iBinnev & Tremain3ll987l) . Thi s is a particular ca se of the distri- 
bution function investigated by ICuddefor3 il99ll) at the limit of 
his parameter — > oo, and thus it is straightforward to derive 
an inversion for the unknown function fE{E), which is scarcely 
an y more difficult th an Eddin gton 's inversion (eq.O (e.g., eg. 24 



oflCuddefordll99lL see also |Deionghelll98dlwrikrnson & Evan3 
ll999tlBaes & Deionghj200a) - 



f(E,L): 



1 



21" 



1 



12/1 (27t)3/2 r(l • 

21" 1 



1 



dE Jo 



(27T)3/2 r(l-Cl')r(l-/3)L2/* 



Jo 



dtf/ 



dif/ d"h 
(E - 1//)" diF 
d"+'/i 



1 d"h 



(E - ij/Y di/^"+' dijj" 



equivalent result. However, strictly speaking, for a ^ I, the inte- 
gral in general diverges and the formula therefore becomes mean- 
ingless although ad-hoc extension of the formal definitions can be 
employed to provide the proper final result. 

For the potential-density pair of equation J4}, we find (G = 
M = a = l) 



47t ^ 



. ^py-2{i-i3)ip^ 



Next, using 



d 

dij/ 



= A- 



(1 - ^"Y 



+ pB 



(1 - 



(1 - l/r")" 
one finds that 

^ =0 

for m < p + 3 - 2/3. Then, equation <AH reduces to 
2'^-' p+1 1 di/' 



(A2) 



nE,L): 



J_ f 



(1 - >P"y 



(27t)5/2 r(l - Q-)r(l -P) L^P Jo (E - if/r di/^"+' 

where j + I = 2(1 - y3)/p > 0. If y > 0, equation fA2l further 
implies that, for < £ < i/' < 1, the integrand is always non- 
negative and so is the distribution function. In fact, it is, at least 
formally, possible to derive the series expression of the distribution 
function for arbitrary p and fS from the integral form, i.e., 



f(E,L): 



p + 1 £"+3/2-/? " Y(,pk + p + 4-2/3) T{k + y) 



(27t)5/2 r(l E^ 



r(pk + p + 5/2- j3) r(r) k\ 



On the other hand, if - 1 < y < 0, equation <A2> also indicates 



that 
d'"h 
dip" 



< 0, 



and subsequently, one can establish that 
lim f{E, L) < 

In other words, there exists a certain value E_ such that f(E, L) < 
for < £ < £ < 1, that is, the corresponding distribution func- 
tion is unphysical. Hence, for the potential-density pair of equa- 
tion ^4), the constant anisotropy distribution function is physical 
only if p < 2(1-/3) and the hypervirial models which have an 
anisotropy parameter /3 = 1 - (p/2) are the models with the maxi- 
mally radially biased velocity dispersions for a given p and a con- 
stant /3. 



(4=0 



.(Al) 



Here, h(ifi) 



'p and /3 = (3/2) - n - a where n is an integer 



and < Q- < 1 . Note that, formally, the choice of an arbitrary non- 
negative integer for n in equation <AH actually does produce the 
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